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We consider the -B — > 7T7T and B — ► 7r/j, 7ru; decays alongside each other, 
taking into account the contributions from all individual penguin amplitudes 
generated by the internal t, c, and u quarks. We argue that three ratios of 
penguin amplitudes, each for a different internal quark, formed by dividing 
the individual penguin amplitude in B — > tttt by the corresponding amplitude 
in B —> up, 7tuj, should be equal. We study the implications of the assumed 
existence of this connection between B — > tttt and B — > irp, ituj. First, ac- 
cepting that in the B — > tttt decays the ratio C/T of the colour-suppressed 
factorization amplitude C to the tree factorization amplitude T is negligible, 
we determine the ratio of individual penguin amplitudes. Then, from the 
B — > irp, ttlo data, we extract the effective (i.e. possibly containing some pen- 
guin terms) tree and the effective colour-suppressed amplitudes relevant for 
these processes, and the corresponding solutions for the factorization ampli- 
tudes. Finally, we argue that the C/T ratio in B — > tttt should be identical 
to its counterpart in B — > irp, ttlo (relevant for pion emission from the decay- 
ing b quark). This constraint permits the determination of C/T and of other 
amplitude ratios directly from the data. Although the \C/T\ ratio extracted 
from the available data still carries a substantial error, it is consistent with 
the expected value of 0.25 - 0.5. 
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1 Introduction 



Decays of B mesons to charmless final states provide us with a lot of information 
concerning the phases of the Cabibbo-Kobayashi-Maskawa matrix. Consequently, 
these decays have been the subject of numerous studies. The problem is to extract 
relevant information from the data, since any such procedure involves serious un- 
certainties resulting from our poor knowledge of the effects of strong interactions. 
Many papers have been devoted to the analysis of B — > PP decays (with P denot- 
ing a pseudoscalar meson), in the hope that the abundance of data will permit the 
determination of several weak and strong parameters involved. 

In recent papers [1] it was shown that the data on B — > niv decays require the 
presence of important nonfactorizable corrections and hadronic interference effects 
if the SM value of 7 ~ 65° and the CP-averaged Bd — > 7r°7r° branching ratios 
recently measured by BaBar and Belle [2] are used. The authors of ref. [1] show in 
a theoretically clean way that the effective tree and colour-suppressed amplitudes 
T and C governing the B — > im decays are roughly equal in absolute magnitudes, 
thus contradicting the expectation that the latter amplitude should be substantially 
suppressed. In order to conform to this expectation, one has to admit that the 
corrections in the effective amplitudes T and C arising from the usually neglected 
penguin contributions are substantial (see also [3] and [4]). Thus, hadronic- level 
effects appear to invalidate the naive expectation of the factorization prescription. 
Indeed, final-state strong interaction effects should contribute to the redefinition 
of the original quark-diagram amplitudes, thus generating the usually neglected 
penguin contributions referred to above (see refs. [5, 6]). 

In general, better extraction of the relevant strong and weak parameters requires 
considering a larger body of data. Recently, several analyses appeared in which 
extraction of the angle 7 of the unitarity triangle (UT) was attempted from a fit 
to all currently available data on B — > PP or B — > PV decays (see e.g. [6, 7, 
8]). Such analyses in the B — > PP and B — > PV sectors are usually performed 
separately from one another. On the other hand, arguments may be given that 
some of the parameters, introduced in these two sectors to take account of the 
effects of strong interactions, are actually related. By combining in a single analysis 
the information from both sectors, one could then hopefully get an additional handle 
on the previously undetermined parameters. 

In this paper, we analyse the decays of B to nn, up, and ircu in an approach mod- 
eled on ref. [1]. The data on these decays should provide sufficient information on 
which to base the analysis and comparison of the size of all factorization amplitudes 
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involved. In principle, there is no need here to use the data on strangeness-changing 
two-body decays of B mesons, a welcome feature since B — > ttK decays exhibit vari- 
ous puzzles (in addition to substantial contributions from electroweak penguins), as 
analysed in [1]. In practice, however, the data on B — > im are still not good enough, 
and we find it necessary to determine the magnitude of penguin amplitude P from 
B + — > 7i + K°. This transition and the related B + — > tt + K*° decay constitute the 
only places wherefrom information from the strangeness-changing sector enters into 
our analysis. The knowledge of \P\ permits the extraction of several B — > ixp^uj 
parameters from the data and the determination of amplitude ratios, such as e.g. 
T/P etc. (i.e. it gives the size of various amplitudes in units of \P\). 

In Section 2, we set out our notation following refs [1, 4] and, using the Summer 
2004 data given by the HFAG [9] we repeat this part of the B — > tttt analysis of ref. 
[4] which is relevant for our purposes. We then present a simple yet illuminating 
formula which expresses the ratio of penguin amplitudes involving loops with differ- 
ent quarks in terms of the ratio of factorization amplitudes C/T and the parameters 
extracted from the data. In Section 3, using the data on B + — > tt + K^° decays, we 
determine the absolute magnitude of penguin amplitudes and formulate our main 
assumption concerning their ratios. In Section 4, we proceed with an analysis of the 
B — > irp and B — > ttuo decays. We find that there are two acceptable sets of solutions 
for the effective colour-suppressed and tree amplitudes. Assuming that the ratio of 
penguin amplitudes involving loops with different quarks is well approximated by 
the formula of Section 2 with C/T = 0, we determine the tree and colour-suppressed 
factorization amplitudes in B — > np, ttuj. In Section 5, we give PV analog of 

the formula given in Section 2, expressing the ratio of penguin amplitudes in terms 
of the ratio of factorization amplitudes in the B — > np, ttuj sector. We combine this 
formula with its counterpart from Section 2 and directly from the data determine 
both C/T and the ratio of penguin amplitudes involving loops with different quarks. 
Our conclusions are contained in Section 5. 

2 Decays B — > tttt 

The B — > tttt amplitudes may be expressed in terms of amplitudes P (penguin), T 
(effective tree), and C (effective colour-suppressed): 



V2A(B 




[T + C] 



(1) 



A(B° d - vr+vr-) 
V2A(B° d - vrV) 



[T + P] 
[C-P], 



(2) 
(3) 
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with: 

P = P C = AX 3 V tc (4) 

f = e ^(T-R b P u ) = AX 3 R b e il (T - V tu ) (5) 

C = e^(C + R b P u ) = AX 3 R b e^(C + V tu ), (6) 

where the rightmost forms involve the definitions of refs [1, 4], from which we omitted 
the contributions due to the exchange amplitudes. In principle, the latter could be 
included by a mere redefinition of P tu . However, the general idea of this paper would 
then require an analogous redefiniton of penguin amplitudes in B — > np, ituj, with 
the relative size of penguin and exchange amplitudes unlikely to be the same as in 
B — > Tin. Only if exchange amplitudes are neglected (as usually done and assumed 
hereafter) and P tu represents just the difference between the top- and up- penguins, 
our approach does not need additional parameters. 

In Eqs (4-5), A = 0.83 ± 0.02, A = 0.224 ± 0.0036, R b = 0.37 ± 0.04, and 7 
parametrize the Cabibbo-Kobayashi-Maskawa (CKM) matrix. T = AX 3 R b T and 
C = AX 3 R b C involve the strong amplitudes of colour-allowed and colour-suppressed 
tree diagrams. Finally, P q = -XfV tq = AX 3 V tq = AX 3 (V t - V q ) with V k de- 
scribing penguin strong amplitudes corresponding to internal k-quark exchanges 
(k G {t, c, u}), i.e. with the full penguin amplitude given by X^V U + X^V C + X^Vt, 
where = V q kV* b , and V is the CKM matrix. For flavour-symmetric final-state in- 
teractions (FSI), the above formulae encompass all elastic and inelastic FSI with the 
exception of those represented by crossed diagrams (see [5, 6]). When the latter are 
taken into account, the amplitudes T and C become mixtures of the colour-allowed 
and colour-suppressed factorization amplitudes [5, 6]. 



2.1 Extraction of hadronic parameters 

Introducing the hadronic parameters of refs [1, 4]: 

de w = -e^P/f = - \P/f\e~ i5 T = -P C /(T - R b P u ) (7) 
xe iA = C/f = \C/f\e^ s c- s f) = (C + R b P u )/(T-R b P u ) (8) 

where Sf, 5^ denote strong phases (in the convention in which the strong phase of 
the penguin amplitude P c is assumed zero), one can derive the following formulae 
relating the parameters just introduced to the branching ratios B and asymmetries 
in B — > 7T7T decays [1]: 

l + 2xcosA + a: 2 

+ - 1 - 2d cos 9 cos -f + d 2 { ' 
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_ d 2 + 2dx cos( A - 9) cos 7 + x 2 

00 " l-2dcos#cos 7 + d 2 1 J 

dir 2rfsingsin7 

" +7r " 1 -2cicos#cos7 + d 2 1 ; 

ix sin(2/3 + 2 7 ) - 2rf cos g sin (2/3 + 7) + d 2 sin(2/3) 

" + "~ ~" 1 -2rfcos^cos7 + rf 2 ' 1 j 



with [9] 



B(B ± -> 7T ± 7r ) 1"o0 

#™ = 2 ^ V^T — = 2-20 ± 0.31 (13) 

+ B{B d -> 7T+7T-) r B + 1 ; 

- 2 ^^"^ = 0.66 ±0.13 (14) 

0(-Dd — > 7T + 7T ) 

= -0.37 ±0.11 (15) 
= ±0.61 ±0.14, (16) 

where, as in original papers [1, 4], the asymmetries are estimated as weighted aver- 
ages, in spite of the BaBar and Belle results still not being fully consistent. Since 
these averages have not changed much in comparison to information available before 
Summer 2004, we believe that performing the analysis of this paper for the BaBar 
and Belle asymmetries separately is unwarranted. 

Assuming f3 = 24° and 7 = 65°, one can determine d and 9 using Eqs (11,12) and 
the experimental values of the asymmetries from Eqs (15,16) [1]. Two solutions for 
(d, 9) are obtained, of which refs [1, 4] accept only the one with small d, the other 
one (with d ss 4.6) being excluded as it leads to complex solutions. 

With the updated values of asymmetry averages, we obtain 

d = 0.52+°; 22 (17) 
9 = +(14112)°. (18) 

In [1] ([4]) the corresponding values were: d = 0.49±ali (O-SlloJo), and 9 = 
+(137l^)° (+(140ll|)°). We determine the errors as in [1], i.e. the errors asso- 
ciated with the specific input parameter (here: asymmetry or an R nn ratio ) are 
estimated by varying its value within la while keeping the other input parameters 
at their central values, with the individual errors thus obtained subsequently added 
in quadrature. 

The solution of Eqs (9,10,13,14) yields 

x = 1.131°;* (19) 
A = -(551^)°, (20) 
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to be compared with x = 1.22±g:ii (i-^-o.ll), and A = -(7lt^)° (-(57l§g)°) in 
refs [1] and [4] respectively. As in [1], we have discarded the second solution for x 
and A (with x ~ 0.96, and A as +33°), since the Acp(B + — > fT°7r + ) asymmetry it 
yields is of the order of 0.2 (using Eq.(33) below), too large when compared with the 
experimental data of [9]: A e ^p{B + — > K°tt + ) = —0.02 ± 0.034. A similar argument 
was originally used in [1] in connection with Acp(B + — > K + tt°). The experimental 
value of the Aqp{B + — > K°tt + ) asymmetry puts a much stronger bound on the error 
in A — 9 (see below) . 



2.2 Ratio of penguin amplitudes 

From Eqs (7,8) one derives: 

P c _ V tc _ _ n C R b de ie 

Since the value of \C/T\ is expected to be of the order of 0.25 only, a good estimate 
of the ratio of penguin terms should be obtained by setting C/T = above: 

# - - (0.1718:") e^s)° (22) 

t u X6 



with 

A-9— (+1641^)° (23) 

where from now on we treat the errors of d, 9, x, and A as independent, adding in 
quadrature the corresponding errors to the moduli and phases of quantities depend- 
ing on these parameters. Since the errors on d, 6, x and A are actually inter-related, 
another possible treatment of errors would be to vary the original four parameters 
(-R+1, -Roo i ^+7i— > ^™+7i—) within their error bars. We have not chosen this more 
involved route since possible errors stemming from the lack of consistency among 
the BaBar and Belle asymmetry results for A^^, ^4™+*- might render its higher 
quality questionable. Error analysis of this paper, with K^l, Rqq , ^+ 7r -, ^™+*- 
varied within their error bars, could be repeated when the input data are under 
better control. 

From Eq. (7) one further obtains : 

T = R, 'T - \ e "* K " (3 ' 6 " ?) e ' <10:i "'' (24) 



and 



1 + V iA ^(1.68±°;^e* (26± - )0 (25) 



R b P u x 
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which demonstrates that the size of the penguin correction term (Rf,P u ) with respect 
to the tree amplitude T is large, as discussed in [1]. 

When small nonzero values of C/T are admitted, the ratio P c j P u receives a 
correction term 

whose inclusion, as shown by the r.h.s. above (obtained by inserting the central 
values of x, d, etc.), increases the errors in our estimate of \P C /P U \ given in Eq. (22) 
by some 20% (for \C/T\ « 0.2). 



3 Size of penguin amplitudes 

Given the current errors in the determination of the \T/P\ ratio (c.f. Eq. (17)), 
further information on the size of amplitudes considered in Section 2 is best obtained 
from the decays of B + -> n + K° (B~ -> ir-R°) and B+ -> ir+K*° (B~ -> ir-R*°) 
as these are expressed in terms of penguin amplitudes alone: 

A(B+^7i + K°) = P' (27) 

A(B + ^n + K*°) = P'p, (28) 

with 

P' = -\l s) V t u-\ { c s) V t c (29) 

P'p = -^Vp.u-X^Vpm, (30) 

where the subscript P ( v ) for B — > amplitudes denotes amplitudes in which 
the spectator quark ends in the final P (V) meson. As already remarked in the 
Introduction this constitutes an implicit use of SU(3), at present necessary, given 
the size of errors in the nn sector. In fact, knowing the absolute size of penguin 
amplitudes well is important in the formulas of Section 4 (see e.g. Eqs (69-72)), with 
the route via Ktt etc. decays yielding much smaller errors. With the size of penguin 
amplitude depending upon the SU(3) assumption in question, the latter affects also 
the extracted values of tree and colour-suppressed amplitudes. 
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3.1 B+ 



7T 



+K° 



Let us consider the B + 
analog of P c ): 



ir + K° and B — > 7r K° decays. We introduce P c ' (the 



P c ' = -Xi s) V tc = -A\ 2 (l - \ 2 /2)V tc . 
Using Eq. (22) for the ratio of Vtu/Ptc one obtains from Eqs (29,31): 



P = E 



l + eR b 



P 



tu 



'tc 



1 -e 



x e 



iA 



d e i6 



(31) 



(32) 



where e = A 2 /(l - A 2 ) = 0.05. 

The CP-asymmetry Acp{B + — > jr+fT ) is approximately 

-2ef sin(A-#)sin 7 



A CP {B + 



(33) 

9 I 7+O.86 
z - 1 1 -0.70) 



1 — 2e^ cos(A — 0) cos 7 

which, together with its experimental value of —0.02 ± 0.034, and x/d 
points towards A — 6 ps (+5lg°)° or (+174 ± 11)°, significantly improving upon 
the value of Eq. (23). Thus, the experimental value of the Acp(B + — > n + K°) 
asymmetry forces A — 9 to be close to 0° or 180°, rejecting the solution with A m 33° 
(for which A — 9 ~ —108°). Consequently, from now on, whenever only A— 9 appears 
instead of A and 9, we shall use the average of the two determinations: 



A- = +(173 ±10)°. 
Thus, the right hand side of Eq. (22) is replaced with 

§ - (0.1718:") e^°. 



(34) 



(35) 



The CP-averaged branching ratio of B + — > n + K° is approximately: 



B(B ± -> n ± K { \K Q )) w |P ( 



./|2 



1 — 2e— cos( A — 0) cos 7 



(36) 



■+ TT^ ^ )) = 

24.1 ± 1.3 (in units of 10~ 6 ), x/d = 2.171^, and A - 9 = +(173 ± 10)°, one finds: 



Using the experimental value of this branching ratio: B exp (B 



IP'I 2 ~ 

Pel ~ 



22.11} 



+1.3 
4 



(37) 



i.e. 



P c ' = -4.7018;^, (38) 

where the strong phase of P' c (originating from P tc ) is assumed zero by SU(3) sym- 
metry with P c . One then finds 

A 



P c = 



1 - A 2 /2 



P' c = -(0.230 ± 0.004)P C ' = 1.08 ± 0.04 



(39) 
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3.2 B+ 



7T 



+K 



*o 



For the description of B + 
P c , we introduce: 

P' 



ir + K*° and B — > 7r decays, in analogy to P c ' and 



-AWp Pitc = -AA 2 (1 - \ 2 /2)V PM . 



(40) 



To proceed we assume that 



Pp. 



\tu 



P,tc 



Ptu 
Ptc 



(41) 



The above equality follows if one accepts that the formation of the final PP or 
PV pair is independent of penguin transition occurring before that formation takes 
place. This should be so if the intermediate su state does not remember how it was 
produced. The relevant penguin-induced decay amplitude becomes then a product 
of a penguin term and the amplitude describing the formation of the final state, with 
the latter amplitude cancelling out in the ratios Vp,tu/Pp,tc and Ptu/Ptc- The above 
assumption is a crucial assumption upon which the rest of this paper is based. 



Assuming Eq. (41), an analog of Eq.(32) follows: 



pi _ pi 

r P — r P,c 



l + eR b 



tu 



tc 



p' 



1-e 



x e 



iA 



d e i6 



3 «7 



(42) 



Consequently, from the experimental branching ratio of B exp (B — > ir K* (K* )) 



9.76i 1 ; 2 2 and using the analog of Eq. (36): 
B(B ± -> n ± K*°(K* )) w \P' P f 



x 



one determines that 



P' 



± 



1 — 2e— cos(A — 9) cos 7 
a 



(43) 



(44) 



Introducing the ratio £ = Pp :C /P c so that the B — > PV amplitudes may be later 
expressed in units of the B — > PP penguin amplitude P c (to be compared with Eqs 
(24,25)), one then finds that 

£ = P PiC /P c = Vp, tc /Vtc = Vp,tu/Vtu = P'p : JP' c = 0.64 ± 0.04 (45) 

with 

p' Pc = £P c ' = -2.99 ± 0.21, (46) 

where we adopted the convention of a vanishing strong phase for the P' Pc penguin 
amplitude (Vpjc)- One then finds 

P P , C = -(0.230 ± 0.004)P; jC = 0.69l°;°i (47) 
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4 Decays B — > Tip, nuj 

Besides the amplitudes Pp(v) (Pp(v),c) already considered in the previous section, 
strangeness-conserving decays of B mesons into a pseudoscalar-vector meson pair 
introduce several further amplitudes: Tp(y) (tree), Cp<y) (colour-suppressed), Sp(y) 
(singlet penguin), etc., of which T P ^ V ), Pp(v),Cp(y) are considered to be the domi- 
nant ones. The Zweig rule suggests that Sp should be negligible. On the other hand, 
Sy does not need to be, in analogy to the situation in the B — > PP sector where 
decays B — > Kr)(r)') seem to indicate the non-negligible size of S' [10, 11]. Since we 
want to restrict our analysis to a group of decays akin to B — > titt, i.e. not involving 
singlet penguin amplitudes, we are left with the following B d ,B + decay channels 
to be considered: n + p~ , n~p + , n°p + , n + p°, n°p°, ty + uj, and ti uj (together with CP 
conjugate processes). By restricting our analysis to these decays, we do not need to 
introduce the additional parameters related to the singlet amplitudes. Our omission 
of the strangeness-changing B — > PV decays is deliberate, since such amplitudes 
seem to exhibit some anomalous behaviour already in the B — > PP sector [12]. 

The respective B — > up, -kuj amplitudes are given by (in sign convention used e.g. 
in ref.[8]) 

A(B° d ^n + p-) = -[f v + P v ] (48) 
A(B° d ^n-p+) = -[fp + Pp] (49) 

A(B + ^n°p + ) = --^[fp + Cy + Pp-Py] (50) 
A(£+-7rV) = -^[fy + Cp-Pp + Py] (51) 
A(B° d ^n°p ) = - l -[C P + Cy-P P -Py] (52) 

A(B + ^n + u) = -^[fy + Cp + Pp + Py] (53) 
v2 

A(B° d ^n°u) = ^[Cp-Cy + Pp + Py] (54) 

where 

P P(y) = P P(nc = A\ 3 V P{vltc (55) 
and Tp(y), Cp(v) involve expressions similar to Eqs (5,6). 

Following the arguments given in [10] and used in previous discussions [11], we 
assume that Py = —Pp, or, more precisely, that: 



V,tc(u) 



-Vp,tc{u)- (56) 
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Then, denning P P)U and P v , u in analogy to Eq. (55) 

P P{v) ,u = AX 3 V P{V)M (57) 

and inserting Vp jtc (u) = ^Ptc{u) from Eq. (45), the amplitudes Tp(y) and Cp(y) may 
be written as 

f v = e^(T v -R b P VtU ) = AX 3 R b e^(T v + CV tu ) (58) 

fp = e i i(Tp-R b P P:U ) = A\ 3 R b e l \T P -iV tu ) (59) 

C v = e^(C v + R b P v , u ) = A\ 3 R b e^(C v -^V tu ) (60) 

Cp = e^(C P + R b P PyU ) = A\ 3 R b e^(C P + £V tu ) (61) 



with Pp u = £P U = —Py :U . In the above equations, the rightmost entries are given 
in the form completely analogous to that used in ref. [1], with Tyrp) = AX 3 R b Ty(p) 
and Cv(p) = A\ 3 R b Cy(p) involving the strong amplitudes Ty(p) and Cy(p) of colour- 
allowed and colour-suppressed tree diagrams. 



Since the amplitude P P = P P)C = —P v is known (Eq. (47)), from Eq. (48) 
and the knowledge of experimental asymmetries and branching ratios for the B% — > 
7r + p~ decay one can determine the magnitude and (relative) phase of amplitude 
Ty. Then, using Eq. (58) with Py u = —Pp tU and the estimate Pp tU = Pp fi P u jP c ~ 
— Pp c a;e^ A_ ^ / (R b d) one can extract the tree amplitude Ty. A similar procedure 
applied to Eqs (49,59) yields tree amplitude Tp. A subsequent use of Eqs (50,51) 
should permit the determination of Cy and Cp. The remaining three equations 
(52,53,54) provide additional constraints/check on the extracted values of colour- 
suppressed amplitudes. We now turn to the extraction of the relevant amplitudes 
from the data. 



4.1 Extraction of tree amplitudes Ty and Tp 

In analogy with Eq. (7), we first introduce the following parameters in the B — > PV 
sector: 

dpe i6p = -^ = -1^1^ = - P ^ p (62) 

1 p — ti b fp u 



e^ P f = - 


\Pp 


T P 


\Tp 


e^ Pv =- 




Ty 


\Tv 



<h-^ = ' — — "'^-- - r • (63) 

ly + ti b F Pu 

For the CP-averaged branching ratio 

B(B° - tt V) = \[B{B» - tt V) + B(B° d - (64) 
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and the CP-asymmetry 

A(B° d - n+p-) - 

one derives: 

B(B° d -> tt V) 
- tt V) 

with (Eq. (47)) 

Solving Eqs (66,67) one gets 
tan Oy - 



B(B° d -> 7r-p+) - B(B° d -> vr+p-) 
B(BS-7r+p-)+B(5g-7r-p+) 



1 + j — COS #y COS 7 



4 . /j . 

^— sin fly sin 7 



1 + i-^ cos ^cos 7 



P P = 0.69 ± 0.05. 



"V 



dy = 



shi7(cos7 ± ^ cos 2 7 + /y — 1 — fcy / sin 2 7) 
sin 6>y sin 7 



v 



where we defined 



A(B° d -> tt V)B(BS - tV) 



4|Pp| 2 



fcy 

With the experimental values: 

B(bS->ttV) = io-Hr.9 



A(B° d -> vr+p-) 



-0.47 



+0.13 
-0.14 



(65) 

(66) 
(67) 

(68) 



(69) 
(70) 

(71) 
(72) 



(73) 
(74) 



taken from [9], after neglecting the correlations with B d — > 7r p + , one obtains the 
following two solutions: 



or 



dy — dy\ 

Qy = Ova 



dy — dy,2 
0y = 6>y,2 



206+ - 025 

U.ZUU_ 022 



= (-34^)° 



(-139lg)°. 



(75) 
(76) 



(77) 
(78) 
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From the CP-averaged branching ratio B(B® — > ti p + ) and asymmetry A(B^ — > 
7r _ p + ) defined in analogy to Eqs (64, 65), using the experimental values 

B(fi d ^7r-p + ) = 13.911? (79) 
A(B° d ^n-p + ) = -0.15 ±0.09 (80) 



from [9], one similarly obtains: 



or 



d P = d PA = 0.174lffi (81) 
6p = e P>1 = (-121D (82) 



dp = d P , 2 = 0.2031^1 (83) 
Op = P , 2 = (-166l|)°. (84) 



Further experimental constraints are given by the parameters S pn and AS pn 
extracted from the time-dependent studies of (P°,P°) — > p ± n T , and providing in- 
formation on the relative phases of effective tree amplitudes Ty and Tp. In terms 
of our amplitudes of Eqs (48,49) and their CP-counterparts, these parameters are 
expressed as follows: 



with 



= (S+- + 5_+)/2 (85) 
AS pn = (S + _-S_ + )/2 (86) 



2 Im A+- 

S+ - = 1 + IA+-P 



(87) 



S -+ = i + |A-+|2 (88) 
' A(P° - P+7T-) rf F e*T - d P e^ 1 J 

A(P° -> p-7T+) ^ dp - dve^v iyUJ 

The four pairs of solutions, i.e. (Pi, Vi), (Pi, V 2 ), (P 2 , Vi), and (P 2 , V 2 ), give the 
four predictions for S p7T and AS pw gathered in Table I. 
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Table 1: Four predictions for S pw and AS p7T 



case (a) case (b) case (c) case (d) 

(PuVJ (P U V 2 ) (P 2 ,Vi) (P2,V 2 ) 

s pn -0.2718:81 -0.0418:81 o.ool™ 4 +0.32+8:8! 

as P7T +o.38l8:i +o.73l8i? -0.7318:18 -o.38l8:i 



Table 2: Relative sizes and phases of effective tree amplitudes 









case (a) 


case (b) 


e 


-ilf v /t 

5 

fp/fy 


c 


(3.nl8:i)e i(+34 -- )0 

(-158l} 7 )° 

i.i8l8:?8 


(2.68l8: 4 8)e t(+139 -? )D 
(-53111)° 

i qy+0.26 
L -° ' -0.22 



Since, according to the data [9]: 

S pw = -0.15 ±0.13 (91) 
ASpn = +0.25 ±0.13 (92) 

cases (c) and (d) may be rejected. Although case (a) is clearly the best, case (b) 
cannot be ruled out: the difference of the two determinations of AS p7T (case (b), Eq. 
(92)) is consistent with zero at (slightly above) 2a. 

In both cases (a) and (b), we have 

e-*f P /P c = (3.6818:*D e*(- 168 -?)°. (93) 

For e-^fv/Pc, the relative phase 5 = kig(f v /f P ) = 6 P - 6 V - 180°, and the 
relative size Tp/Ty = dyjdp of tree amplitudes, one finds the results gathered in 
Table 2. These solutions should be compared with 8 ~ —22° and \Tp/Ty\ ~ 1.46 
obtained in the favored fit of ref. [8], and corresponding to our case (b). 

One of the essential differences with ref. [8] is the fact that in our approach the 
effective tree amplitudes T P and T v do not correspond to the tree amplitudes T P 
and Ty of the factorization picture. Instead, the effective amplitudes T P and T v 
involve substantial corrections to the factorization terms, due to the presence of the 
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Pp,u (Pv,u ) part of the penguin amplitude (Eqs (58,59)). Using Eqs (62,63) one can 
determine the tree amplitudes: 

£ - h(^-^) (94) 



P, ' 'I', dp 

Z> p-V 

u ° 

P c \ P c dy 



From Eq.(94), using the estimate (see Eq.(22)) 



P c d 

and Eq.(34), for both cases (a) and (b) one obtains: 



R ^ « - (96) 



^ = (2.40^)e J (- 157 -3)°. (97) 

Similarly, from Eq.(95) one gets: 

case (a) ^ = (2.25^;^)e J ( +58 ^)° (98) 

case(b) ^ = (3.91^;^)e l ( +150 -i")° (99) 

For case (a) one finds |7V/7>| = 0.94l^, while for case (b): \T V /T P \ = 1.63±8;||. 
If the B —> p and B — > 7r formfactors are similar, one expects (see [8]) that the ratio 
of |7V/Tp| should be approximately equal to the ratio of f n / f p ~ 0.63, as Ty (Tp) 
involves a weak current producing ^ (p 1 * 1 ). Thus, case (a) seems favoured again. 

For both cases (a) and (b), however, one also estimates from Eqs (24,97) that 
\T P /T\ « 0.671q;25, which disagrees with the simple expectation (c.f. [8]) of \T P /T\ « 
f p /f w = 1.59 (while \ f P /f\ = id/dp = 1.91±g;||). Still, one has to keep in mind 
that the above estimates are based on Eq. (96) which neglects terms of order C/T. 



4.2 Extraction of colour-suppressed amplitudes Cp and Cy 

In analogy with Eqs (62,63), we introduce the following parameters involving colour- 
suppressed amplitudes C P and C v : 

ype tFp = -Tf^rTT^ (100) 

y v e^ = -^e- (101) 

J-P + Oy 
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The CP-averaged branching ratio for the B + — > 7T + p° decay and the correspond- 
ing asymmetry are given by 



2 



B(B + ^n + p°) = (4 + ^ + 40087^^) — (102) 

V Vp Vp J 2 

A(B + ^n + p°) = -4sin 7 ^^ = ^. (103) 

1 ' VP 4 + 4 r + 4cos 7 £2 ^ 



Vp 'VP 



Using the experimental numbers of 



B(B + ^n + p°) = 9.1 + 1.3 (104) 
A(B + ->ttV) = -0.19 + 0.11 (105) 



one finds two solutions: 
sol. (PI) 



y P = y P , 1 = 0.195±g;g|l (106) 

r = r P) i ee (+23 + 14)° (107) 



and 

sol. (P2) 



y P = y P , 2 ee 0.149l°;°£ (108) 
r P = r Pi2 ee (+16311?)°. (109) 



From Eq.(100) one has: 



3 n 



(110) 



Putting the estimates of yp, aV etc. into Eq.(llO), one obtains the values of 
Q-E-e' 11 given in Table 3. 

For the £> + — > 7r°p + decays, one obtains formulas completely analogous to (102) 
and (103), with yp — > yy, and Tp — > Ty. Using the experimental branching ratio 
and asymmetry: 

£(£+ ^ 7r°p + ) = 12.0 + 2.0 (111) 
A(S + ->7r°/o + ) = +0.16 + 0.13 (112) 
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Table 3: Effective color-suppressed amplitudes Cp from B + — > 7r + p° decays 



e-*C P /P c 


case 


(a) 


case 


(b) 




sol. (PI) 
sol. (P2) 


(5.6218:21) 
(i.6i^:l) 


e i(-175±\l)° 
e i(-17tS)° 


(i.n±S:S) 
(6.1518:2?: 


) e* 


-5t 9 9 )° 



Table 4: Effective color-suppressed amplitudes Cy from £> + — > n°p + decays 



sol. (VI) 7.5318$ e l ( +17 -">r 
sol. (V2) 2.47lJi? e i ( +117 --)° 



one finds two solutions: 
sol. (Vf) 



and 

sol. (V2) 



y v = y v , 1 = 0.16518:821 (H3) 
T V = T VA = (-211H) (114) 



y y = y y , 2 = 0.13018:811 (115) 
Tv = T v ,2 = (-163111)°. (H6) 



From Eq.(lOl) one has: 




(117) 



Putting the estimates of yv, dp etc. into Eq. (117), one obtains the values of 
%^e-^ given in Table 4. 

By comparing the experimental branching ratio for B® — > 7r°p° and the bound 
on — > -k°lu (from [9]) with the predictions of all combinations of entries in Tables 
3 and 4, one finds that only cases (a,P2,V2) and (b,Pl,V2) may be admitted. We 
shall refer to them as Solutions I and II respectively. The corresponding predictions 
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Table 5: Branching ratios for — > 7r°p° and B° d — > tt°uj decays 



cxp 



Sol. I 



Sol. II 




5.0 ±1.8 1.68 



+0.83 
-0.51 



0.86 



+0.90 
'-0.49 



< 1.2 



2.28 



+0.81 
-0.65 



2.31 



+0.87 
-0.80 



for the branching ratios of — > 7r°p° and B® — > tt°c<j are compared with the data 
in Table 5. 

Discrepancies with experiment observed in Table 5 suggest that the assumptions 
(in particular the SU(3) assumption of Section 3 and/or possibly Eq. (41)), which 
lead to the value of Pp given in Eq. (68) thereby affecting the extracted size of 
colour-suppressed B — > PV effective amplitudes, might not be wholy adequate. 
One needs here a way of estimating the size of Pp in the strangeness-preserving 
sector, which would be both sufficiently precise and less assumption-dependent. 

Solution II, with yp m 0.195, is fully consistent with the information gained 
from the branching ratio B(B + — > tt + uj) = 5.9 ± 0.8, which yields yp = 0.20 ± 0.02. 
Solution I, with yp m 0.149, agrees with yp determined from B + — > tt + cj at 2a. 
Thus, the B + — > h + uj branching ratio favors Solution II over Solution I. We recall 
that it is just the opposite case with the values of S pn and AS p7T which favor case 
(a) (hence Solution I) over case (b) (Solution II) by 2a. Since for the B® — > 7i°p° 
branching ratio, as Table 5 shows, the difference between experiment and theory 
is 1.7a (2.1a) for Solution I (II), one concludes that Solution I describes the data 
slightly better than Solution II. 

The colour-suppressed factorization amplitudes may be estimated from 




(118) 



(119) 



and 



Pc 




(120) 



y v dp 




,*(A-(9) 



(121) 
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Table 6: Color-suppressed factorization amplitudes Cp and Cy obtained for C/T = 



Sol. I Sol. II 


c P 


Re(C P /P c ) = 0.1518:1 C P /P C = 2.40li:°Je*(- 173 -) 
Im(C>/P c ) = -0.311};°} 


Cy 


C v /P c = 2.05l};ie i ( +82 -")° 



The values of Cp and Cy obtained for both Solutions I and II assuming C/T = 
(Eqs (119,121)) are gathered in Table 6. 

Interestingly, with the central value of \T P /P C \ being 2.40, Solution I is consistent 
with a small value of C P /T P) while for Solution II the C P /T P ratio is of the order of 
I. On the other hand, given the central T v /P c value of 2.25 (3.91) for Solution I(II) 
respectively (Eqs (98,99)), it is Solution II for which Cy /Ty seems to be smaller. 
One has to remember, however, that in our calculations we used the values of x and 
d determined from the averages of not fully consistent asymmetries in the 7T7t sector. 
Furthermore, our estimate of errors in the determination of P c /P u did not include 
the errors due to nonvanishing C/T. As remarked earlier, for \C/T\ m 0.2 these 
corrections may increase the error of P c /P u by 20%, affecting the ensuing discussion 
correspondingly (see Section 5). 

5 Extraction of C/T and P c /P u 

In the analysis performed so far, the ratio C/T of the factorization amplitudes in 
B — > 7T7T decays has been assumed negligible. It turns out, however, that one can 
actually determine the value of C/T directly from the data, provided one is willing 
to make an additional very plausible assumption. Namely, we observe that the 
amplitudes C, T in B — > nn decay and the amplitudes Cy, Ty in B — > up, ircu 
transitions are due to the same process, namely a decay of b quark into a pion and 
a light quark. The difference between the two colour-suppressed amplitudes C and 
Cy (and between the two tree amplitudes T and Ty) should be due only to the fact 
that the amplitude for the recombination of the freshly produced light quark with 
the spectator quark depends on whether the two recombine into a pseudoscalar 
(in B — > 7T7r) or a vector meson (in B — > 7ip,7iuj). However, this dependence on 
the recombination amplitude should cancel in the ratios, i.e. in C/T and Cy/Ty. 
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Consequently, we may assume that 

% = % (122) 

T Ty y ' 

We now recall Eq. (21), which correlates the ratio P c /P u with the size of C/T. We 
seek a similar connection for the B — > np, nu sector. To this end, we observe that 
using the expressions (59,60) for the effective amplitudes in Eq. (101) we can write: 

Now, Eqs (62,63) may be rewritten as 

TL = f:( Rb ^-JL e -»A (124) 



I = « + ^ * ) ' (125) 



From Eqs (123,124) we determine 



P c V p c Vv dp I 



Dividing Eq. (126) by Eq. (125) we obtain 



r Rh E* + J- e - iT v + -1- 

Cv_ nb P c + yv C + dp 



e -i6 P 



T v ' J_ e -iO v _ Rh Eu 

dy^ Ub P c 



(127) 



The above equation may be rewritten in the form completely analogous to Eq. (21), 
namely: 

dp Ty dy 

In the denominator above, the first two terms partially cancel. By assuming that 
Cy/Ty is so small that the third term may be neglected, we obtain a counterpart 
of the previous estimate of P c j P u given in Eq. (22): 

F = - i P -»,*i,-uv K O-^SSK^, (129) 

dp y v 

which, despite the approximation involved, is consistent with Eq. (22), and thus 
with a large value of P u as compared with P c (and a small value of C/T). 

If C/T is assumed equal to Cy/Ty, Eqs (21,128) may be solved for C/T with 
the final result: 

_L e -t0 P i -L e - ir v _ £ e *(A-0) 
- - ^ ^ d . (130) 



T _Lp-i8v _ l p -it 

d v d e 



20 



Table 7: Extracted values of C/T 



Solution I Solution II 

C/T (0.4718$) e*( +47 -")° (0.75181) e^ 30±2 ^ 



Let us take the central values of the parameters and discuss the denominator first. 
Solutions I and II differ in their values for the parameters of the pair (d v ,9 v ). For 
Solution I (II), the first term in the denominator has an absolute value of around 4.9 
(4.2). The second term has an absolute value of around 1.9. The sum of the terms in 
the denominator, with phases taken into account, has an absolute value of 6.8 (4.3) 
for Solution I (II) respectively. As for the numerator, our previous considerations 
uniquely determined the values of the three numerator terms. In particular, the 
absolute value of the first term is equal to 1/dp = l/<ip,i ~ 5.7, that of the second 
term is 1/yv = l/yv,2 ~ 7.7, while for the third term it is equal to x/d pa 2.2. It is 
therefore non-trivial that with the central values of phases taken into account, the 
sum of these terms is not large and has the absolute value of around 3.2, leading 
to the central value of \C/T\ for Solution I being 0.47. The sum of the three terms 
in the numerator is most sensitive to the value of angle IV. If IV is set at its la 
deviation value of —177°, the absolute value of the numerator becomes 1.3 only. For 
Solution I, the value of \C/T\ would then become equal to 0.2. When all of the 
errors are calculated, one obtains the values given in Table 7. 

The determinations of Cy/Ty — C/T given in Table 7 may be compared with 
the central value of C v /T v = 0.91e i24 °(0.52 e - i68 °) for Solution I (II) obtained in 
Section 4 for C/T = 0. Thus, the previously obtained central value of Cy/Ty 
gets significantly reduced (increased) for Solution I (II). Although in Solution II the 
central value of \C/T\ is now quite large it is also compatible with \C/T\ of order 
0.25. Better data are clearly required. 

Other estimates of C/T also lead to values of order 0.5. For example, in ref. 
[1] arguments in favor of C/T = 0.5 x e t290 ° are given. Similarly, in their recent 
SU(3)-symmetric fit to all B — > PP decays, Chiang et al. [3] obtain the value 

\C/T\ = 0.461°$. 

With the central values of \C/T\ in Table 7 significantly larger than the expected 
value of around 0.25, the original estimate of P c /P u , obtained in Eq. (22) upon 
assuming C/T = 0, could be substantially affected. Solving Eqs (21,128) for P c / P u , 
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Table 8: Extracted values of C P /T P 



Solution I 



Solution II 



Re{C P /Tp) 
lm(C P /Tp) 
Cp/Tp 



-n 24+ - 23 
-0.031 036 



(0.2518 



+0.31^ 
21; 



-0.35 



+0.40181 
+0.00+8:1 



(0.4018 



0.40\ 
26 J 



,<(0tS)° 



one obtains 

p 1 — - e~ i61 

^ = R b d e ld ^ ™— , (131) 

where 

k ~ (132) 
Numerically, for Solution I one finds: 

= (0.2118:S) e^3)°, (133) 



,+21 n 



which still bears resemblance to (0.1718:81) e <(16 ^)° f Eq. (22). 
On the other hand, for Solution II one obtains: 

# = (0-71±S:S) e l(+20 - )O , (134) 

with error estimates so large that they admit small values for both \P C /P U \ and 
\P U /P C \. Again, there is here a strong dependence on T v , with larger values of 
\P C /P U \ attained when T v is set at its la deviation value of —148°. 

From Eqs (94,118,131), one can further determine the corresponding values of 
Cp/Tp. They are gathered in Table 8. 

From Tables 7 and 8 we see that it is Solution I which prefers smaller central 
values of both Cy/Ty = C/T and Cp/Tp. With present errors, however, both 
Solutions I and II are still compatible with \C/T\ and \C P /T P \ of around 0.25. 



For completeness, we have also calculated the ratio \T V /T P \ obtaining 0.9618: 
(0.8818:26) f° r Solution I (II) respectively (to be compared with the value of f w / f p 
0.63 expected in [8]). 
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6 Summary 



In this paper we performed a joint analysis of the B — > txtx and B — > 7rp, 7tc<j decays 
with the aim of studying the effects of the presence of two independent superpositions 
of penguin amplitudes on the possible values of colour-suppressed and tree factor- 
ization amplitudes. Our analysis assumes that the formation of the final PP or PV 
pair is independent of the penguin transition occurring before that formation takes 
place. This constitutes a crucial assumption of our approach. The analysis yields 
two sets of solutions for the effective colour-suppressed (Cy,Cp) and tree (Tp,Ty) 
amplitudes in the B — > np, nu transitions, with one solution weakly favoured over 
the other one. 

Assuming the C/T ratio in B — > tttt to be negligible, we estimated the ratio of 
the two superpositions of penguin amplitudes, using it subsequently to determine 
the values of Cy,T v ,Cp,T P from the data. This procedure yielded two sets of 
numerical estimates for Cp/T P and Cy/Ty. 

By imposing the condition of equality for the ratios of C/T and CyjTy we de- 
termined the value of C/T directly from the data. The two solutions obtained are 
compatible both with a value of \C/T\ of around 0.25 and with the estimates from 
literature yielding \C/T\ ~ 0.5, with errors still of the order of 0.3 — 0.4. The 
corresponding solutions for P c /P u and Cp/Tp have been given as well. One of the 
solutions is preferred as it yields smaller central values of both C/T = Cy/Ty and 
Cp/Tp. Discrimination between the solutions, and a more precise determination of 
C/T, require better data. When such data become available, a well-defined value 
may be extracted for C/T along the lines similar to those presented here and com- 
pared with expectations and other estimates, providing us with more information on 
the C/T ratio and the expected connection between penguin amplitudes in B — > nn 
and B — > np, ircu decays. 

This work was supported in part by the Polish State Committee for Scientific 
Research (KBN) as a research project over the period of 2003-2006 (grant 2 P03B 
046 25). 
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